Random Antiferromagnetic SU{N) Spin Chains 
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We analyze random isotropic antiferromagnetic SU (N) spin ciiains using tiie real space renormalization 
group. We find that they are governed at low energies by a universal infinite randomness fixed point different 
from the one of random spin-| chains. We determine analytically the important exponents: the energy-length 
scale relation is f2 ~ exp (— L"^), where tp = 1/A'^, and the mean correlation function is given by dj ~ 
(—1)'"-' / \i ~ j\ , where <j> = 4/N. Our analysis shows that the infinite-A'^ limit is unable to capture the 
behavior obtained at any finite TV. 
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The identification of several universality classes of disor- 
dered quantum systems in one dimension (ID) has seen enor- 
mous progress in recent years. Prominent among those is 
the case of random antiferromagnetic spin-i chains, which 
have been shown to be controlled by an infinite random- 
ness fixed point (IRFP) at low energies fl, 2V Many prop- 
erties of this so-called random singlet phase have been ob- 
tained, e. g., the spin susceptibility x ^ ^/T\og T, and 
the spin correlation function Cij — (S^ • Sj) is such that its 
mean value Cij ^ ( — 1)* ■^ / \i — j\ , while the typical one 
CijL ~ exp(— -^/ji — j|). Further studies of ID spin sys- 
tems have uncovered a wide variety of behaviors such as var- 
ious Griffiths phases L2.^^^i^^ and large spin phases 

It is the purpose of this paper to extend these analyses by en- 
larging the symmetry group from SU (2) to SU (N). We have 
several motivations for this. The inclusion of orbital degrees 
of freedom often leads to the enlargement SU (2) -^ SU (N). 
The strong spin-orbit interaction in rare-earth elements locks 
spin and orbital moments into a large multiplet with degen- 
eracy N, whose description requires the enlargement from 
SUi2) to SUiN) 1 12]. Recently, a realization of a self- 
conjugate SU{4) spin chain has been proposed in a pillar ar- 
ray of semiconducting quantum dots, where the symmetry- 
breaking effect of the intra-dot electron-electron interaction is 
minimized due to the peculiarities of the dot potential 1 13]. 
Several other possible realizations of enlarged symmetry have 
been considered in the literature (see, for example, Ref. 11411 '). 
In any of these cases, the effects of disorder are clearly of in- 
terest. Furthermore, the large- A^ limit of SU{N) spin models 
is of considerable interest. In this limit, many models can be 
solved by saddle-point methods and 1/A^ corrections can be 
obtained in a controlled manner ] 15]. The hope behind this 
approach is that the physics of A^ = 2 is at least qualitatively 
captured as A^ — > cxd. The solution of a sequence of models 
as a function of N, though rarely possible, can determine the 
validity of this idea. We will show that the random antiferro- 
magnetic SU (N) chain provides just such a solution. Inter- 
estingly, the low-energy physics at finite A^ is never captured 
at infinite A^. We will show, however, that for any finite A^, 
the system is governed by a new universal IRFP with charac- 



teristic exponents which we calculate and depend only on the 
group rank. 

We will focus on the following Hamiltonian 



H 



J2J^r^ 



i+lj 



(1) 



i=l 



where Ji are positive independent random variables dis- 
tributed according to Po{J), and the components of Ti, Ff 
(a = 1, . . . , N'^ — 1), are the generators of a representa- 
tion of SU{N). We will confine our analysis to totally an- 
tisymmetric irreducible representations. These correspond to 
Young tableaux with one column and Qi lines ] 16]. They are 
conveniently expressed with the help of auxiliary fermionic 
operators q„ (a = 1 ,..., A^) through Ff = cl^T'^^c^^, 
where F°^ are the generator matrices of the fundamental 
representation of SU{N). The fermions obey the constraint 
Eq=i ^ia^ia ~ Qi- ^^ considered the cases where the Q^s 
are random and Qi — Q — const. 

To treat the Hamiltonian Q, we generalize the real-space 
renormalization group method introduced by Ma, Dasgupta 
and Hu ] 1 ] . Our generalization is reminiscent of the treat- 
ment of random ferro- and antiferromagnetic spin chains ] 1(1]. 
We first find the largest bond energy of the system, fi — 
maxjAi}. We define A^ as the energy difference between 
the ground and first excited multiplets of the z-th bond. As 
Ji > 0, it can be shown that the ground state multiplet is 
represented by a vertical Young tableau with Q lines, where 
Q = Qi + Q,+i,ifQ, + Q,+i < Af,andQ = Qi + Qi+i-N, 
if Qi + Qi+i > N. The energies of ground and excited multi- 
plets can be calculated from the Casimir's with the usual trick 

2r, • r,+i = (r, + r,+i)' - rf - rf^^. The vaiu^of the 

Casimir's of the relevant tableaux is given in Ref. ] 17]. We 
then decimate that bond by keeping only the ground multiplet 
and renormalizing the neighboring interactions in the follow- 
ing fashion. If Qi + Qi+i = N, the bond ground state is a 
singlet and is thus removed from the system. The new effec- 
tive coupling between the neighboring spins F^-i and ri+2 
is, by second-order perturbation theory. 



J 






(2) 



Otherwise, the spin pair Ti and r^+i is replaced by a new ef- 
fective spin r, which belongs to a totally antisymmetric rep- 
resentation with Q lines as given above. It connects to the 
spins Fi-i and ri+2 through renormalized couplings given in 
first-order perturbation theory by 



Ji~i = CiJi~i, and Ji+i = (1 - ^i) J, 



J+I7 



(3) 



respectively, where £_i = Qi/ {Qi + Qi+i), if Qi + Qi+i < 
N, and ^,; =^ {N - Qi) / {2N - Q, - Qj+i), otherwise. We 
point out the similarity with the case of the random chain with 
both ferro- and antiferromagnetic interactions LIOJ . where 
both 1st and 2nd order decimations are generated. Unlike the 
latter, however, here the active (i.e., not yet decimated) spin 
clusters are always vertical tableaux and the procedure always 
maintains a totally antisymmetric spin chain. Moreover, the 
renormalized couplings are always smaller than the original 
ones. Thus, at every decimation step, the energy scale fi is 
lowered. 

An important feature of the decimation procedure is that it 
does not privilege any specific representation. When a spin 
pair (Qi, Qi+i) is decimated out, the new effective spin is 
never equal to any of the original ones (Q j^ Qi, Q ^ Qi+i). 
Thus, after an initial transient, each one of the A^ — 1 totally 
antisymmetric representations is equally probable, even if we 
start with Qi = Q = const, (except for some special fine- 
tuned cases dealt with later). We have confirmed this numer- 
ically, as will be shown later. We have also checked that the 
distribution of representations becomes uncorrected with the 
distribution of couplings. We thus focus on the flow of the 
coupling distribution, P (J) = P {J; fi), as the highest scale 
il is decreased lll|2|]. We take = 1 initially. 

As will be shown later, similarly to the random spin-i 
chain, P (J) always flows to an extremely broad distribution. 
We are thus justified in neglecting the numerical prefactors 
in Eqs. (|2j and (|3}, which are always less than unity and ir- 
relevant asymptotically |2]. Furthermore, there are a total of 
{N — 1) ^ possible decimation processes, all of them equally 
likely. Of these, A^ — 1 are 2nd order, each with probability 
p = 1/ {N — 1), and the others are 1st order, with probability 
q = 1 — p. Thus, we can write a flow equation in the useful 
logarithmic variables F = — In 17, and C, = In (il/ J) Ii2j| 



— p(C; r) = — p(C; r) + qp{C; r)p° + pp°p ® p, (4) 

where p(C; r)dC = P (J; n) d J, p° = p(0; T), and p p = 
/ dCidC3p(Ci; r)p(C3; r)6 (C - Ci - Cs). The first term on 
the right-hand side is due to the fact that C, changes when F 
increases. The second one, absent in the random spin-i chain, 
is due to 1st order decimation steps and only ensures the nor- 
malization of p. The last one is due to 2nd order steps, which 
strongly renormalize p broadening it. 

If Po{J) is not extremely singular, the flow Eq. © has only 



one stable fixed point solution ^ 



P*{J;n) 



pT 

a /n 



9{n-j), 



(5a) 
(5b) 



with a = l/(pF) = - (iV - 1) / In Q. The fixed point distri- 
bution Q broadens indefinitely in the limit fi ^ 0, rendering 
the renormalization procedure increasingly more precise, and 
asymptotically exact |2]. The system is thus governed by an 
IRFP. 

The relation between energy and length scales can be de- 
termined by finding the fraction of active spin clusters nr 
at the energy scale F ||2|]. If F is increased by dF, a frac- 
tion driY = {2p + q) nrp{0; F)dF of active spin clusters is 
decimated. Thus, close to the fixed point, where p(0; F) « 
P*(0;F)ta 



.1/^ _ 



= [l^(l/f7)]l/^ 



(6) 



where tp = p/{p + l) = 1/N. This type of 'activated' dynam- 
ical scaling, corresponding to a dynamical exponent z ^ oo, 
arises here with an unexpected exponent ip. When N — 2, 
we recover the usual form found in the random spin-i chains 

Si. 

In order to check the validity of the approximations leading 
up to Eq. 0, we have numerically implemented the full pro- 
cedure. The data were generated by decimating chains with 
lengths up to 10^, averaging over 100 realizations of disor- 
der All the initial spins belong to the fundamental represen- 
tation (Qi = 1, \/i). We analyzed several initial distributions 
Po{J). Chains A, B, and C had uniform distributions in the 
interval a; < J < 1, with x — 0.9, 0.5, and 0, respectively. 
In chains D, E, and F we used initial power-law distributions 
Po ^ J"'^, with (3 = 0.3, 0.6, and 0.9, respectively. In 
Fig. n we show the fraction of first order decimation steps 
as a function of the energy scale fl, for the symmetry groups 
SU{3) and 5/7(4). As anticipated, it tends asymptotically to 
q = {N — 2) / {N — 1). The figure also shows the fraction 
of Q = 1 spins in the two cases, and the fraction of self- 
conjugate (Q = 2) spins in the SU{4) chain. They all tend 
asymptotically to 1/ {N — 1), as expected. 

In Fig. |2j we plot nr as a function of F for SU{3) and 
SU{A). By fitting the asymptotic behavior, we confirm the 
universality of the exponent t/) = 1/A^, as predicted by Eq. (|6|l. 
We point out that tp converges in a logarithmic manner, thus 
a more precise determination of tp demands the decimation 
of longer chains than the ones studied here. We see that as 
N increases, so does the number of decimations needed for 
a given decrease in energy scale. This 'delayed scaling' can 
be understood by realizing that only 2nd order processes are 
effective in lowering the energy scale, and these become less 
frequent as N increases. 

There are other IRFP's in addition to the one analyzed 
above. For example, the self-conjugate SU{2k) spin chain 
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Figure 1 : The fractions of first order decimation steps, of spins in the 
fundamental {Q = 1), and in the self-conjugate (Q = 2, only for 
S(/(4)) representations, all as a function of fi. For clarity, we only 
show data for chains A (solid lines) and E (dashed lines) (see text). 
The filled (open) symbols refer to the ^[/(S) (5'(7(4)) chains. The 
data error is about the size of the symbol. 



o 




Figure 2: Fraction of active spins nr as a function of the energy scale 
F, for the SU{?,) (main) and S'(/(4) (inset) chains A to F. The data 
error is smaller than the symbol size. 



(with integer fc > 1) flows towards an IRFP with t/" — 1/2, 
since Qi = fc, Vi, and only 2nd order decimation steps oc- 
cur, hke in the random spin-i chain. Although these chains 
are gapful 1 18], they are unstable against the introduction of 
weak disorder, due to the topological nature of their ground 
state, as explained for the random J\ — J2 Heisenberg chain 
in Ref. ^. More importantly, this V' = 1/2 IRFP is unstable 
against the introduction of Q ^ fc spins. For a small concen- 
tration rii of such spins, the system will initially be governed 
by the V' = 1/2 IRFP, until the energy scale T ^ n^ ' is 
reached. Below that scale, the renormalization flow veers to- 
wards the IRFP of Eq. Q, with the characteristic exponent 
-0 = l/2fc. Similar IRFP's exist for other SU {N) chains, but 
they are equally unstable with respect to the introduction of 
'defect' spins. 



We can easily calculate the asymptotic behavior of thermo- 
dynamic quantities using Eq. (|6j |2]. Since P{J) becomes 
very broad at low energies, the active spins are approximately 
free at a low temperature T — Vt, whereas the decimated ones 
do not contribute, since they are frozen in singlet states with 
excitation energies much greater than T. Hence, we find that 
the entropy density a 



cific heat c 



InT) 






(-InT) ^'"'^ and the spe- 
Furthermore, it can be eas- 



ily shown that the magnetic susceptibility of a single SU (N) 
spin is Curie-like, from which it follows that for the whole 
system x (T) - nr/T = 1/[T (- InT)^/'^]. 

We can also obtain the asymptotic behavior of the spin cor- 
relation function Cy = (F^ -Tj) |2|. Spins belonging to the 
same cluster develop O (1) correlations; otherwise, they are 
weakly correlated. Therefore, such spins dominate the mean 
correlation function. To find Cij , we need the probability that 
any two well-separated spins Fi and Tj are rigidly locked 
in the same spin cluster T when |z — i| ^ n^^. First, we 
need to find V {t;nr), the probability to find a spin cluster 
r composed of t original spins at scale nr- After many dec- 
imations, any spatial correlations between Q's and J's have 
vanished and any remaining bond is equally likely to be dec- 
imated. The fraction of clusters with t spins at scale np is 
firV {t; nr). When dNdec decimations are performed, nr de- 
creases by dnr = — (2p + q) dNdec and 



d[nrV{t)] = 



-2Vit)+qY,rih)Pit2)St, 



+ t2:f 



dN, 



dec 7 



where the two terms on the right-hand side give the fraction of 
decimated and added clusters with t spins and we suppressed 
the nr dependence of V to Ughten the notation. Hence, 



nr- 



dV (t) 1-p 



dnr 



l+p 



vit)-J2'P{h)r{t~ti) 



(7) 



whose solution is V{t:nr) ~ npexp(— inp) in the limit 
nr -^ 0, with 7 = (1 - p) /(I + p) = 1 - 2/N. Finafly, 
the probability that Ti and Tj are active in the same cluster is 

equal to (tnr) ~ ( np"'' j , yielding. 



Oi' 



(-1) 



»-j 



(8) 



with 4> = 4:/N. The typical correlation function, however, is 
very different. Following |^, we note that it involves many 
factors of ,/ decimated at various scales e"'^. The scaling 

behavior is dominated by the smallest factor O ( e^'^l'^-'l 

where k ~ 0(1) yielding typical correlations 



\C, 



^3 Hyp 



exp 



{-k\^-Jf) 



(9) 



Fig. |3] shows numerical results for the mean correlation 
function a~ for groups SU{2), S'[/(3) and S'[/(4), aver- 
aged over 200 realizations of disorder for chain lengths up to 



S£/W, (|)=1.01 

su(3),^=^.3^ 

Sf/('2), (|)=1 .99 




above [20]. In particular, the energy-length scale exponent is 
the same. In that case, N is the number of phases meeting at 
the multicritical point. 

In conclusion, we have identified in random antiferromag- 
netic SU (N) chains an infinite randomness fixed point with 
exponents different than the ones previously found in spin- 
i chains. An important question which we leave for future 
study is the stability of this phase against the introduction of 
anisotropy. 
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Figure 3: Mean correlation function for groups SU{2), SU{3) and 
S'f/(4). The power-law dependence is evident. The exponents <jf> — 
1.99 ± 0.03, 1.31 ± 0.03, and 1.01 ± 0.03 are obtained by fitting 
the region r > 10, and are in excellent agreement with the analytical 
value (j) = 4/N (see Eq. <H)). 



L = 10^ and open boundary conditions. The numerical proce- 
dure consists in completely decimating a chain, and counting 
the fraction of spin pairs that become strongly correlated at the 
distance |i — j| [7]. Excellent agreement with the analytical 
prediction of is obtained. No significant dependence on ini- 
tial disorder strength was observed, confirming the universal 
behavior. 

In the large-iV limit, the mean correlation function decays 
extremely slowly. In this limit, the fraction of 2nd order pro- 
cesses is very small and the mean number of spins in a cluster 
diverges at low energies, all of them being strongly correlated. 
A 1/iV expansion of Eq. (|8} leads to |Cy | ^ 1/ \n{\i — j\). 
Incidentally, this is the same behavior observed numerically 
in random ferro- and antiferromagnetic spin chains [19.1 . This 
is no surprise, since both systems are dominated by similar 
1st order decimations whose clustering rules are the same as 
N ^ oo. Therefore, they are both described by Eq. Q with 
p — 0, hence the logarithmic dependence of the mean corre- 
lation function. As far as we know, this analytical explanation 
has not appeared before. However, we should stress that the 
asymptotic region governed by the IRPP is reached at energy 
scales which decrease with the increase of N, since the 2nd 
order processes become increasingly rare. Therefore, in the 
infinite-iV limit the universal behavior described above dis- 
appears and a direct infinite-A^ approach fails to capture the 
physics at any finite N. 

Interestingly, some multicritical points of random antifer- 
romagnetic spin S chains have been shown to exhibit a struc- 
ture that is very similar to the generic SU{N) IRFP described 
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